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A METHOD OF COMPUTING THE TRANSIENT TEMPERATURE OF THICK WALLS FROM 
ARBITRARY VARIATION OF ADIABATIC-WALL TEMPERATURE AND 

HEAT-TRANSFER COEFFICIENT 1 

By P. R. Hill 


SUMMARY 

A method of calculating the temperature of thick walls has 
been developed in which are used relatively new concepts, such 
as the time series and the response to a unit triangle variation 
of surface temperature, together with essentially standard 
formulas for transient temperature and heat how into thick 
walls. The method can be used without knowledge of the 
mathematical tools of its development. The method is partic- 
ularly suitable for determining the wall temperature in one - 
dimensional thermal problems in aeronautics where there is a 
continuous variation of the heat-transfer coefficient and 
adiahatic-wall temperature. The method also offers a con- 
venient means for solving the inverse problem of determining 
the heat-flow history when temperature history is known. 

A series of diversified problems were solved ,by exact analysis 
as well as by the new method. A comparison of the results 
shoios the new method to be accurate. The labor involved is 
very mndcst in consideration of the nature of the thick-wall 
temperature problem. Limiting solutions for the “infinitely 
thick” wall and for walls so thin that thermal lag can be 
neglected were also obtained. 

INTRODUCTION 

In aeronautical applications, external surfaces are heated 
by the impact and friction of the air. For cases in which 
the structural temperatures never reach equilibrium, the 
transient temperatures of the surfaces often govern the 
design; and it is necessary to be able to predict these tem- 
peratures. 

Literature on transient temperatures in thick walls dates 
from the classical works of Fourier. Perhaps the most ex- 
tensive work on the subject is given in reference 1. Most 
literature giving the solution to the transient temperatures 
in thick walls is based on the premise that the temperature 
history of one or more principal surfaces is known or given. 
Only a limited amount of literature is available relative to 
transient temperatures in thick walls under the influence of 
forced convection. The forced-convection equation for 
heat transfer in aeronautical applications is g=h(T aa —T), 
which states that the rate of heating g is proportional to 
the difference between the adiabatic-wall temperature T aw 


and the wall temperature T. The coefficient of propor- 
tionality is the heat-transfer coefficient h. In the classical 
problem of the convection heating of a thick wall, h has 
been assumed to be constant. In the usual aeronautical 
application, the fact that h varies with time is the source of 
the difficulty in obtaining a solution. 

The thick-wall case treated in this paper is the one gov- 
erned by Fourier’s classical partial differential equation 

bT j &T 

bt k bA 

In the case governed by this equation, the wall is composed 
of a homogeneous material and the temperature gradients 
and heat flow parallel to the surface are negligible. In one 
boundary relation for this case, the convective heat rate is 
equated to the heat absorbed by the wall or to the product 
of the conductivity and the temperature gradient in the wall 
at the heated surface; that is, 

Since h occurs as a product with T in this boundary equation, 
the usual procedures of operational calculus do not apply. 
When solutions for the temperatures of thick walls have 
been necessary in aeronautical work, the method generally 
used has been to divide the thick wall into a number of 
slabs in order to make a step-by-step numerical integration 
of Fourier’s equation of heat flow. Since steps in both 
distance and time must be taken, the procedure is tedious 
and time consuming unless the use of a high-speed auto- 
matic computing machine is resorted to. If it is necessary 
to do the work without the use of such equipment, a method 
introduced by Schmidt (ref. 2) wherein some of the calcu- 
lations are accomplished graphically may be used to reduce 
the labor to some extent. This method is known as the 
Schmidt plot method. 

In the present paper a simple method is developed for the 
calculation of the temperature history of the surfaces of a 
thick wall or of any plane within the wall. The procedure 
is to select from a table a set of coefficients which depend 
on the physical properties of the wall. These coefficients 


i Supersedes NAOA Technical Note 4105 by P. R. HU1, 1957. 
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and other data are substituted into explicit algebraic for- 
mulas to determine the temperature of the heated wall 
surface. If the heat-transfer coefficients are known, no 
guess or iteration procedure is required. As can be seen by 
the results of the example problems presented, the accuracy 
can be as good as is desired. For equal time-step sizes, the 
method is more accurate than more laborious numerical 
methods. 

The simplicity of the results depends on two factors: 
One is the suppression of the variable x representing the 
distance into the wall by using an integrated form of Fourier’s 
equation and assigning a value of x corresponding to the 
heated surface. The other is a mathematical device known 
as the time series introduced by Tustin (ref. 3). The time 
series is defined in appendix A. Reference 3 also introduced 
various manipulations of the series. The multiplication of 
two series is an important manipulation by means of which 
specific results can be generalized. Other writers (ref. 4, 
for example) have also presented various manipulations of 
the series. 

The present paper is divided into two parts, analysis and 
application. The section on analysis includes a treatment 
of the determination of the temperature history for the 
special cases of the thermally thin wall and the infinitely 
thick wall as well as for the wall of intermediate thickness. 
The inverse problem of determining the heat flow correspond- 
ing to a known temperature history is also discussed. Al- 
though the method was set up for the purpose of predicting 
wall temperatures in engineering applications, it has also 
been found to be suitable for research applications wherein 
the transient skin temperature is measured and the heat-flow 
and heat-transfer-coefficient histories are deduced. Appen- 
dix A gives background material pertaining to the use of time 
series that may be an aid to a study of the analysis. Ap- 
pendix B gives a summary of analytical temperature and 
heat-flow formulas used either as a basis of analysis or used 
in the solution of examples to test the accuracy of the present 
method. In the section on application the computing 
formulas are reviewed and several examples of their use are 
given. Because of the explicit nature of the temperature 
formulas, it is not necessary to study the analysis to use the 
results. 

SYMBOLS 

A dimensionless coefficient to determine heated- 

surface temperature history 

b slope of adiabatic-wall temperature with respect to 

time, °F/hr 

B dimensionless coefficient to determine inside-wall 

temperature history 
c specific heat, Btu/lb-°F 

e base of natural logarithms 

F{x,t) any function corresponding to reference slope 

y=(lfd)t 

F v (t) any function corresponding to arbitrary control 

line y=y(t ) 

F&(t) any function corresponding to unit triangle control 
line 

heat capacity of wall, pci, Btu/(sq ft) (°F) 


Tjh heat-transfer coefficient at x=Q, Btu/(hr)(sq ft) 

(°F) 

h heat-transfer coefficient at x=l, Btu/(hr)(sq ft) 

(°F) 

H heat-transfer number, A 3^/1 6 67 

heat-transfer number for infinitely thick wall, 
Zh-yJnk8/8K 

k diffusivity, Klcp, sq ft/hr 

K conductivity, (Btu)(ft)/(hr)(sq ft)(°F) 

1 wall thickness, ft 

M memory coefficients, dimensionless 

to term designating time in multiples of basic inter- 

val 5 

n term number in infinite series 

2 instantaneous heat-transfer rate due to uniform 

temperature rise of heated wall surface of 1° 
in time 5, Btu/(hr)(sq ft) 

q m average heat-transfer rate from time (to— 1)5 to 

to 5 due to uniform temperature rise of heated 
wall surface of 1° in time 5, Btu/(hr)(sq ft) 

2 a heat-transfer rate corresponding to unit triangle 

variation of surface temperature 
2 a, m average heat-transfer rate from time (to— 1)5 to 

to 5 due to unit triangle reference temperature 
variation of heated-wall surface, Btu/(hr)(sq ft) 
2 ,(t) heat-flow history at heated surface due to unit 

temperature step of heated surface 
r radiation rate, Btu/(hr)(sq ft) 

R radiation term, r Sir 2 / 160, °F 

R a radiation term for infinitely thick wall, 3R-y,lwkS/8K, 
°F 

t time, hr 

T heated-wall-surface' temperature, °F 

T s value of step in wall-surface temperature, °F 

T aa adiabatic-wall temperature or effective boundary- 

layer temperature, °F 

T t temperature of inside (unheated) surface or of any 

. plane within wall, °F 

V(x,t) temperature response to unit step in T aa , °F 
x distance through wall, ft 

y ordinate of control line or altitude of triangle 

a„ positive roots of auxiliary equation in analytical 

solution of wall temperature 
5 basic time interval in time series 

t] ratio of heat-transfer coefficient at cooler wall 

surface to heat-transfer coefficient at heated 
wall surface 

0 difference in temperature between heated surface 

and any other plane due to unit triangle varia- 

tion of heated surface, °F 

0 T difference in temperature between heated surface 

and any other plane due to uniform reference- 
temperature rise of 1° in time 5, °F 
p weight density, lb/cu ft 

r dummy time variable, hr 

Subscript or superscript: 
to index denoting term number 

A dot over a symbol denotes the derivative with respect 
to time. 
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ANALYSIS 


TEMPERATURES ON OUTSIDE SURFACE 


The wall considered in this paper is composed of a homo- 
geneous material, and the temperature gradients and heat 
flow parallel to the surface are negligible. The transient 
temperatures of the heated or outside surface of the wall 
are determined by means of Fourier's equation which governs 
the heat flow through the wall : 


bT, &T 


( 1 ) 


The wall properties expressed by the diffusivity k are con- 
stant. The boundary conditions are given by the state- 
ments that the flow of heat at the unheated face of the wall 
(where x is taken as zero) is zero, that is, 



and that the heat transferred to the heated face is given by 
the relation 

h(T aw ~T)=K(bT/bx) s „ t (3) 

For convenience, the initial temperature is taken as zero at 
zero time: 

T=0 (t= 0) (4) 

Since both h and T are functions of t and occur as a product 
in equation (3), the transform procedures of operational 
calculus do not apply. The problem may be stated in 
another form by means of an integral equation expressing 
the heat balance at the heated surface. 

Let q,(t) be the heat-flow history at the heated wall surface 
at x=l due to a unit step in that wall-surface temperature 
at f=0. Then the heat flow q(t) due to a temperature 
history T(t) of the heated surface that is continuous and is 
zero when £=0, but is otherwise an arbitrary variation, can 
be expressed by Duhamel’s integral in the form indicated in 

the following equation, f represents ^jT{t) and r is a dummy 

time variable: 

2 (*)=J o ‘ T(t-T)a,(r)dr (5) 

A heat balance is formed at x=l by equating the integral in 
equation (5) to the left member of the boundaiy condition 
expressed by equation (3) : 

h(T aa -T)=j‘ T(t-r)q,(r)dr (6) 

The determination of T from equation (6) establishes the 
solution. 

The wall temperature T is determined, in general, from 
equation (6) for thermally thin, thick, and infinitely thick 
walls. The method first discussed is based on thick walls, 
and modifications of this method are introduced for the 
special cases of thermally thin and infinitely thick walls. 

SOLUTION FOE THERMALLY THICK WALLS 

Time series. — In the calculation of wall temperature T 

G2GG97 — 00 08 


for thick walls, the right-hand member of equation (6) is re- 
placed by the product of two time series. As explained in 
detail in appendix A, a time series is the value (here the 
Crdinate) of a function of time at successive equal incre- 
ments of time 5. Thus, any variation of wall surface tem- 
perature can be expressed as the series 

T=T U To, T 3 ; . . . T m (7) 

When a wall surface has a unit triangle variation of surface 
temperature, the surface temperature increases from 0° at a 
constant rate to a value of 1° at the time f=5 and decreases 
at a constant rate to the value 0° at t= 25. The heat flow 
corresponding to a unit triangle variation of surface tem- 
perature can be expressed by the series 

2A=gA.l, ffA.2, gA,3, • • • ffA.m (8) 

The product of equations (7) and (8) gives the instan- 
taneous values of heat flow due to the temperature varia- 
tion given by equation (7) and could be used to eliminate 
the integral in equation (6). However, a slight variation 
or refinement of the method is introduced which has been 
found to increase the accuracy of the results without in- 
creasing the labor involved. 

If to represents the term number of a time series, the prod- 
uct mS represents the corresponding time. The refinement 
consists in averaging the heat flow from the time (to— 1)5 to 
the time mS. Let q&, m represent the average heat flow over 
this interval due to a triangular variation in surface tempera- 
ture. Then the heat-flow history due to the triangular 
temperature variation can be represented by the series 

2A=2A.l, gA.2, 2A.3, • • • 2 a, m (9) 

where q&, i is the average for tho time 0 to 5, 2 a, 2 is the average 
for time 5 to 25, and so forth. 

The heat-flow history expressed as the average heat flow 
over successive increments 5, but due to the temperature 
variation (7), is given by. the product of equations (7) 
and (9): 

2= (2\, Ti, T%, . . . Tm) (2a, 1 , 2a, 2 , 2a. 3 j • • • 2A,m) (10) 

Such a multiplication actually gives the result by forming 
the proper superpositions, as demonstrated in appendix A. 

In order to adjust the boundary condition expressed by 
equation (3) or the left-hand side of equation (6) to represent 
an average flow of heat q m over the interval 5, the average 
flow of heat from the boundary layer is approximated by 
the mean of the values at the beginning and end of the in- 
terval. Thus, for the interval ending at t=mS, has the 
value 

2<n=2 lAxiCFaio — TOm+Aw-iCTatt — T)m-l] (11) 

If radiation is important to the problem, it may be in- 
cluded. Let r m be the rate of heat radiation per unit area 
at the time toS. Equation (11) may then be written 

\h m (.T aw — T)„— r m +h a . 1 (T ao - T)„_i— r„_J (11a) 
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With or without radiation, the heating history, or varia- 
tion of may be expressed by the series 

2=2n h, h, ■ • • 2» (12) 

wherein each term has the value given by equation (11) or 
(11a). Equating the average heat flow given by equation 
(12) to that given by equation (10) yields 


2d 2*. 23, • • • Zm= (T u T a , T 3) . . . q A , 2 , 2 a, 3, . . . 2a, m ) 

( 13 ) 

In order to evaluate equation (13), the product in the 
right-hand member is expanded by algebraic multiplication 
and values of both members of the equation pertaining to 
equal time are equated. The following set of equations 
results (for simplicity, the radiation terms are not written in): 


— [Aj ( Taw — T) i-\-k-o ( Taw — T)o] — TigA.l 

2 $2 (Ta * — T)i-\-hi ( T aw — T) i] = T a q A , i+ 2) 2 a, a 


(14) 


2 [h-m(,T aw — T) n +^ m _i(jP as , — T) m _i] — T r ,2A,)+T M _i2A.a+ • • • +2\2 a.„ 
Equations (14) can be rearranged to obtain the equivalent equations for T: 

rp hjTgiC'i j~ ApTa;.. Q 

1 22A.1+AI 


rp hji q t<". a 1 hi 


Tgto t l — hiTi 
22a. 1+^2 


— 2Ti2A,a 




( 15 ) 


rp (hTq n ) m -±XhT aw —hT)w- 1 —2(q A ' 3 T m -i+(i A jT m - a -\- . 
“ 22a, i AK, 


• • +2a,wTi) 


The values of q A , m must be derived. 

Average heat flow due to unit triangle variation of surface 
temperature. — In order to obtain the average heat flow due 
to a unit triangle variation of surface temperature, the 
average heat-flow rate q m due to a uniform increase in the 
sin-face temperature of 1° in each unit time 5 must be 
obtained. The average is taken over the time 5 by inte- 
grating the instantaneous heat-flow rate from t=(m — 1)5 to 
t=m8 and dividing by 6. This determination is carried out 
in appendix B. The result is 


- (A . . . \ 


( 16 ) 


In equation (16), A m and A„_i are the summations 




v* k& 

-ra( 2 n— L) 2 __ 




A m _i=S 


kS 

___ 


" =1 jf(2a-l)< 


(17) 


and 6 is the heat capacity of the wall per square foot per °E 
and is the product of weight density, specific heat, and wall 
thickness: 

6— pci (18) 

The average heat flow due to a unit triangle variation of 
surface temperature q A , m is obtained by the superposition of 


the heat flows q due to throe linear variations of wall tem- 
perature as follows (for further details, see the development 
of equation (Al) in appendix A) : 


2 A.m=2»“ 22„- 1 +2m-a (1°) 

Expanding equation (19) by substituting for q m from equa- 
tion (16) gives 

( 20 ) 

Substituting in equation (20) successive values of m, starting 
with m=l, ignoring any parenthetical group in which nega- 
tive subscripts appear, and collecting terms result in the 
following equations : 

8(?/V . . . \ 

^-AS\8 +Al ~ A ) 

2A,2=^5^“-g +A 2 — 3 Ai+2A°^ 

f (21) 


_ 8 Q 

Q&. m=== . H ** 3-4jn— i“f* 3A m - 2 -<4 W — 3) 
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For convenience, the quantities in parentheses can be tabu- 
lated. The quantities in parentheses usually retain sig- 
nificant values after the completion of the temperature 
triangle which created them. For this reason, in accordance 
with the notation of reference 4, they are called memory 
terms and are designated by the symbol M. With this 
notation, equations (21) become 

- 8<?,, "I 


established : 

M l =£+A l -A 0 

Afa = — — 3A1+2A4 
A/ 3 — A3 — 3-4j-f-3a4i — Ac 


( 22 ) 


- 


■M-m — 3 - 4 m _i + 2 3 J 


^■ z ~A5 Mz 


- - 8G 

^ 9r M-m 


(21a) 


Hence, the following equations for memory terms are 


Obviously, the values of M are combinations of the values 
of A. Values of A and M were computed for a range of 
values of kS/P, and the memory coefficients M are listed in 
table I. The value of M decreases with increasing term 
number and sooner or later further terms can be neglected. 

Resulting temperature formulas. — Equations (21a) give 
the value of ga, m sought to complete equations (15); there- 
fore, equations (21a) are substituted into equations (15). 


The result can be simplified by dividing through by 1 66/A8 and letting 


TT 

(23) 

1 \ 

If radiation is important, the appropriate terms are included by using equation (11a) rather than equation (11). Since 
equations (15) are being divided through by 16G/AS, the radiation term R is defined as 



(24) 


With the substitutions of equations (23) and (24) in equations (15), the final results, including terms for radiation, are 


rp HlT aa 'l-\-HQT av! ' 0 — R\ — f? 0 

MH-H: 

rp H 2 T aw , 2 + (HT aw — HT) 1 — M 2 Ti — R 2 — Ri 

Mi+H 2 

rp &iT aW 'Z + (HT aw HT) 2 — M 2 T 3 — M 3 T 1 — R 3 — R ~2 

M.+H3 


' 1 




(25) 


T • . ■ -M m Ti~R n -R> 

m T 7 i tt 


M 1 +B m 


SOLUTION FOE INFINITELY THICK WALLS 

General considerations. — If a wall is thermally very thick 
and is heated rapidly so that the unh eated side experiences 
little heating, it is convenient and accurate to assume that 
the wall is infinitely thick. The same formulas, equations 
(26), are used to compute the wall surface temperature. 
However, instead of the values of M for a particular wall 
or diffusion number, the values of M which are used are 
always a fixed set of numbers which are now derived. The 
values of H and R are also changed. 


Determination of average heat flow due to unit triangle 
variation of surfaoe temperature. — The determination of the 
heat flow due to a unit triangle variation of surface temper- 
ature of an infinitely thick wall depends upon the instanta- 
neous heat flow into the surface due to a unit rise in surface 
temperature in unit time. From page 110 of reference 5, 
the instantaneous heat flow is equivalent to 22&/7 ffirk. 
Since the heat-flow rate is proportional to the surface- 
temperature slope, the instantaneous heat transfer due to 
unit rise of surface temperature in the time 5 is 2 K-yfij^wk- 
This expression is integrated with respect to t between the 
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TABLE I.— VALUES OF M AND B 


\ tell* 

0.01 

0.02 

0. 05 

0. 1 - 

0.2 

0. 5 

1. 0 

2.0 

6. 0 


M 

1 

0.09280547 

0.13124675 

0.20751933 

0.29347652 

0.41495535 

0.64728198 

0.85683719 

1.02954117 

1.15146420 

2 

-.01592290 

— .02251838 

-.03560496 

-.05051887 

-.07836507 

-.22817155 

-.51142924 

-.82682809 

-1.06920821 

3 

-.02683275 

— .03794724 

-.06002546 

-.08722820 

-.14861561 

-.30052454 

-.31662013 

-.20127699 

-.08224604 

4 

-.00957529 

— .01354153 

— .02166582 

-.03696546 

-.07338883 

-.08405219 

-.02634648 

-.00142616 

-.00000036 

6 

-.00553399 

-.00782683 

-.01320029 

— .02628367 

-.04463344 

-.02447704 

-.00223431 

-.00001026 

-.00000000 

0 

-.00373390 

-.00528487 

-.00983261 

-.02025561 

— .02724662 

-.00712803 

-.00018948 

-.00000007 


7 

-.00274116 

-.00389296 

-.00809228 

-.01579610' 

-.01663399 

-.00207577 

-.00001607 

-.00000000 


8 

-.00212390 

-.00304409 

-.00695594 

-.01233887 

-.01015501 

-.00060449 

-.00000136 



9 

-.00170890 

-.00249270 

-.00608368 

-.00964054 

-.00619962 

- .00017604 

-.00000012 



10 

-.00141412 

-.00211864 

-.00535620 

-.00753254 

-.00378485 

-.00005126 

-.00000001 



11 

-.00119015 

-.00185533 

-.00472749 

-.00588550 

-.00231065 

— .00001493 

-.00000000 



12 

-.00103005 

-.00166301 

-.00417648 

-.00459860 

-.00141064 

-.00000435 




13 

-.00090058 

-.00151702 

-.00369098 

-.00359308 

— .00086120 

-. 00000 L27 




14 

-.00079788 

-.00140175 

-.00326234 

-.00280743 

-.00052576 

-.00000037 




15 

-.00071529 

-.00130717 

-.00288362 

-.00219357 

-.00032097 

-.00000011 




16 

-.00064810 

-.00122687 

-.00254891 

-.00171393 

— .00019595 

-.00000003 




17 

— .00059290 

-.00115666 

-.00225306 

-.00133917 

-.00011963 

-.00000001 




18 

-.00054712 

-.00109383 

-.00199166 

-.00104635 

-.00007303 

-.00000000 




19 

-.00050883 

-.00103658 

-.00176041 

-.00081756 

-.00004459 





20 

-.00047651 

— .00098373 

— .00155609 

- .00063880 

-.00002722 





21 

-.00044899 

-.00093448 

-.00137549 

-.00049912 

-.00001662 





22 

— .00042532 

-.00088827 

-.00121584 

-.00038998 

-.00001015 





23 

-.00040477 

— .00084473 

-.00107473 

— .00030471 

-.00000619 





24 

-.00038676 

t- .00080356 

-.00094999 

-.00023808 

-.00000378 





25 

-.00037081 

-.00076455 

— .00083973 

-.00018603 

-.00000231 





26 

-.00035655 

-.00072753 

-.00074227 

— .00014535 

- .00000141 





27 

-.00034368 

-.00069237 

— .00065612 

-.00011357 

- .00000086 





28 

-.00033196 

-.00065895 

-.00057997 

-.00008874 

- . 00000053 


, 



29 

— .00032119 

-.00062716 

-.00051265 

-.00006933 

-.00000032 




, 

30 

-.00031121 

-.00059693 

-.00045315 

-.00005417 

-.00000020 






0 

1 

1.00000000 

0.99999992 

0.99956262 

0.98873183 

0.92596579 

0.69945453 

0.45623855 

0.24814441 

0.00099056 

2 

-.00000015 

-.00019241 

-.02166158 

- . 12553206 

-.31355030 

-.48643195 

-.41618829 

-.24630216 

-.09999000 

3 

-.00002494 

-.00333067 

— .06963676 

- . 17636953 

- . 23787104 

-.15098748 

-.03666380 

-.00182900 

-.00000046 

4 

-.00033480 

-.01195026 

-.09012904 

— .14882949 

- . 14587796 

-.04396967 

-.00310842 

-.00001315 

-.00000000 

5 

-.00144190 

— .02167828 

-.08988006 

-.11749067 

-.08906596 

— .01280454 

-.00026361 

-.00000009 


6 

-.00344276 

-.02904775 

-.08284990 

— .09193126 

-.05437468 

— .00372885 

-.00002236 

-.00000000 


7 

— .00595602 

-.03361673 

-.07435643 

-.07184409 

- .03319564 

— .00108689 

-.00000190 



8 

-.00854573 

— .03599626 

-.06609621 

-.06613647 

- . 02026587 

-.00031622 

-.00000016 



9 

-.01092537 

-.03687746 

-.05854671 

— .04386205 

— .01237227 

-.00009209 

-.00000001 



10 

-.01296010 

-.03679326 

— .05179171 

-.03427135 

— . 00755325 

-.00002682 

- . 00000000 



11 

— .01461453 

-.03611078 

-.04579381 

-.02677771 

-.00461124 

-.00000781 




12 

-.01590635 

-.03507183 

-.04048318 

— .02092260 

- .00281516 

— .00000227 




13 

-.01687671 

-.03383230 

-.03578599 

-.01634774 

-.00171865 

-.00000066 




14 

-.01757370 

-.03249148 

-.03163302 

-.01277320 

-.00104923 

-.00000019 




15 

-.01804423 

-.03111189 

-.02796174 

-.00998026 

-.00064055 

-.00000006 




16 

-.01833036 

-.02973243 

-.02471645 

-.00779802 

-.00039106 

-.00000002 




17 

-.01846823 

-.02837687 

-.02184779 

-.00609293 

-.00023874 

-.00000000 




18 

-.01848809 

-.02705931 

— .01931206 

-.00476067 

-.00014575 





19 

-.01841484 

-.02578769 

-.01707064 

-.00371972 

-.00008898 





20 

-.01826874 

— .02456607 

-.01508936 

-.00290638 

-.00005432 





21 

— .01806620 

-.02339609 

-.01333804 

-.00227088 

-.00003316 





22 

-.01782042 

-.02227782 

-.01178998 

-.00177434 

- .00002025 





23 

-.01754197 

— .02121043 

-.01042159 

-.00138637 

-.00001236 





24 

-.01723931 

-.02019254 

-.00921203 

-.00108323 

— .00000755 





25 

-.01691918 

-.01922244 

-.00814285 

-.00084638 

-.00000461 





26 

-.01658695 

— .01829828 

-.00719776 

-.00066131 

-.00000281 





27 

-.01624686 

-.01741811 

-.00636236 

-.00051671 

— .00000172 





28 

-.01690229 

— .01658000 

-.00562392 

-.00040373 

- .00000105 





29 

-.01555588 

— .01578203 

-.00497119 

-.00031545 

-.00000064 





30 

-.01620972 

-.01502236 

-.00439422 

-.00024648 

— .00000039 






* For additional value of kt/P, see table n. 


limits (m— 1)5 and mo. Dividing by 5 gives, for the 
average heat-flow rate over the interval 5 terminating 
at mS, 



[to* 4 — (to— l)* 4 ] 


(26) 


slope function to that of the unit triangle input function is 
accomplished by substituting equation (26) into equation 
(19): 

4 jr 

■== [to* 4 — 3(to— 1)«+3(to— 2)«— (to— 3) w ] 

3-M<5 


The usual superposition required to change the result of the 


(27) 
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Resulting temperature formulas. — If the bracketed quan- 
tity in equation (27) is designated as the memory coefficient 
M m , then 

3(m-l)«+3(m— 2)«— (m— 3)«] (28) 

A dimensionless heat-transfer coefficient (suggested by eq. 

(27)) is defined as 
\ 

(29) 


the following heat balance: 
hm(Ta U —T) m —r„+h m - i{T aa —T) m -!— 

r m-i=jKOT)„-(OT) m -j] (33) 

Solving for T m gives 

^i?aMU+^)?au>,o — h 0 T 0 -{-T 6 0 T 0 — ri — r 0 
T, = - 


and a coiTesponding term for radiation is defined as 


R , 


3r-^irkS 

~8K~ 


(30) 


T>= 


lkT a u.2-\rhlTaa,l — G X T X — — T x 




(34) 


The substitution of equation (27) into equations (15) again 
results in equations (25) . Hence, equations (25) are used to 
obtain the heated-surface temperature of the infinitely thick 
wall as well as of walls of intermediate thickness, except 
that M, H, and 2? for infinitely thick walls are defined by 
equations (28) to (30). Inspection of equations (28) to 
(30) indicates that the wall material properties and time-step 
size are expressed by equations (29) and (30), while the mem- 
ory terms are invariant with wall properties or step size. 
Substituting successive integers for m from 1 to 20 into 
equation (28) gives the following corresponding values of M: 


Mi =1.0 

M 2 = -0.171573 
M 3 = -0.289 129 
M«= -0.103176 
M b = -0.059630 
M # = -0.040234 
M 7 = -0.029536 


M s = -0.022885 
M fl =— 0.018412 
M 10 = —0.015232 
M n = —0.012874 
M u = -0.011069 
M i3 = -0.009650 
Mu— —0.008611 


M,8= -0.007580' 
M u = -0.006807 
M„= —0.006157 
M,8=— 0.005605 >- 
M 19 =— 0.005130 
M m = -0.004719 


These values of M, along with values of and 2?„ from 
equations (29) and (30), can be used in all problems wherein 
the wall is so thick relative to the heating rates and times 
involved that the wall behaves as though it were infinitely 
tliick. 


THIN WALLS 

When a wall is thermally thin, the temperature drop 
through the wall becomes negligible and the problem is 
simplified by assuming that all interior temperatures are 
equal to the surface temperature. The heat absorbed by 
the wall during any time interval S must be equal to the 
gain in enthalpy or total heat during this time. Hence, 

Q m -iT m -i) (32) 



Any variation of 6 with temperature is accounted for by 
equations (34). If the wall properties do not change over 
the temperature range covered, obviously, O is a constant. 
If O is considered to be constant, equations (34) can be 
derived from equations (25) as follows: As the diffusion 
number kS/l 2 becomes large, all values of A approach the 
value A/8, and from equations (22) it is seen that the only 
memory terms not identically equal to zero are Mi and M 2 , 
which have the values A/8 and — A/8. Eliminating M from 
equations (25) and utilizing definitions (23) and (24) yield 
equations (34). 

In equations (34) the terms h 0 , T aw , 0 , and To have been 
retained since in the thin-wall problem, unlike the thick-wall 
problem, it is convenient for T 0 to have any value. These 
equations have considerable advantage because the need 
for temperature extrapolation is reduced, if not eliminated. 
The equations tend to give accurate results and, as is shown 
subsequently, are suitable for the use of relatively large time 
increments. 

INSIDE TEMPERATURES 


If the heated wall surface is called the outside surface, 
the temperatures at other parallel planes may be called 
inside temperatures. In particular, this paper is concerned 
with the inside surface temperature. According fo the 
notation used in this paper, the inside surface is designated 
by xfl= 0, the outside surface by xfl— 1, and other planes by 
values of x/l between 0 and 1. 

Consider a wall, initially at zero temperature, which has 
the heated surface x=l raised at a reference temperature 
slope T—{1/S)t, while the surface x=0 is insulated. The 
difference in temperature 0 T between the heated surface and 
any plane x is shown by equation (BID of appendix B to be 


# '=?s{.^fc=9 cos [ <2 "- I) 5f]- 

{<2»-i)§g 


«° ( i)* _i 

S(2a-l) 3C ° S| 


_ (2n _i, 5 ^ 


*} 


\ 


Equating the average rate of gain in enthalpy as given by 
equation (32) to the average rate of heat transfer through 
the boundary layer as given by equation (11a) results in 


(35) 
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A set of terms is defined, to represent the summations in 
equation (35) : 


„ 16 Z s ^ (-l)”" 1 

B °~^ (2n — l) 3 003 


[(2»-D |f] 


„ 16 p ^ (-l)*- 1 

Bx ~~T? ks ^ i (2a— l) 3 C0S 


— (2n-l)>££i 


-2(2n-l)>^a 


4Ii 


(36a) 


_ 16 P ^ (— l)— 1 

*8§ (2a-l) 3 C0S 


[< 2n - ulf] 




variation of the heated surface would be given by the time 
series d t =8 lt 0 2 , 0 3 , . . . 8 m . In addition, any general tem- 
perature of the heated surface can be represented as 1 = T h 
T 2 , T z , . . . T m . The temperature difference due to this 
general temperature variation is obtained by formally 
multiplying these two time series, and the inside t emperature 
is obtained by subtracting the product from the heated- 
surface temperature T. Hence 


Tu 2 = — (8iT 2 +d 2 Ti) 

Ti,z~ Tz — (®t ?3 -j- 02 Tj + ^3 Tl) 


y (40) 


T i , a =T„-(8 1 T a +8 i T n - 1 + . . . +8 m 7\) J 


If the temperature difference across the entire wall is sought, 
then x/l—0 and equations (36a) become 


t> 16 l 1 ^ 

-°o — ~s Ti 2-t 


T? L V"\ 


B.r 


t> £2 __ 'V' 

7 .x 


(_l)n-l 

P 


(2a— l) 3 

2k5 


(— I)"" 1 . 

— (2n— 

1)5 TH 3 

(2a- 1) 3 

4 P 

(-D- 1 . 

(2a — l) 3 e 

-2(2 n- 


(-l)- 1 . 

— m(2n 

■jP Jk 
-1>* ~ 75 

(2a— l) 3 C 


4 P 


(36b) 


Regardless of whether equations (36a) or (36b) are used, 
equation (35) may be written 


d T — B 0 — B m 


(37). 


By the usual superposition, in order to change the result 
of the reference slope function to that of a unit triangular 
input function, the temperature difference between the plane 
of the heated wall and the plane being considered is 


8 — Bo — B ra — 2(Z?o — 7? m _i) + Z? 0 — B m -z (38) 


Computation of T tm does not depend on prior computation 
of T u j, T u a, and so forth. 

If the temperature distribution through the wall is required 
in a form which is analytical with respect to x and can be 
differentiated or integrated analytically with respect to x, 
the procedures outlined in appendix: C should be followed. 

CALCULATION OF HEAT FLOW FROM TEMPERATURE 
HISTORY 

TEMPERATURE HISTORY OP OUTSIDE SURFACE KNOWN 

If the heated-surface temperature history is known, the 
heat-flow history can be determined by substituting equa- 
tions (21a) into equation (10) as follows: 

on 

(Tl> T2} Ti> . . . T m ) (Mi, M s , M s , . . . M m ) (41) 

Multiplication shows that the mth term is given by 
on 

2„=g (M 1 T m +M i T„- 1 + . . . +M m T 1 ) (42) 

The heat-flow history can be determined readily from a given 
temperature history of the heated surface by moans of 
equation (42). The average heat flow over a small interval 8 
can be assumed to give the instantaneous rate at the center 
of the interval. 


Assigning integral values to m and simplifying give the 
equations corresponding to successive values of time: 

0i= — Z?i~l-Z?o 

0j= — B 2 -p 213 1 — Bo 

8 3 =-B 3 +2B 3 -B l „ (39) 

0i — — R i4-2R 3 — B 2 _ 


&m=— B m -\-2B m -i — B m -2 J 

A suflicient number of values of 0 for practical purposes are 
given in table I. According to equations (39), the tem- 
perature difference between the plane of the heated surface 
and any other plane due to a unit triangle temperature 


TEMPERATURE HISTORY OF INSIDE SURFACE KNOWN 

If the temperature history of the outer surface or of a 
plane near the outer surface is known, the feasibility of 
accurately determining the heat flow is excellent. If the 
temp era tine history of the inside surface for a thermally 
thin wall (Jct/P large and hl/K small) is known, it is also 
feasible to determine the history of the heat flow into the 
outer surface. However, if the wall is thermally thick, 
relatively small changes in temperatures at the inside surface 
may make it difficult to reconstruct the temperature history 
and heat flow at the outer surface. 

A rearrangement of equations (40) may be used to deter- 
mine the outside-surface-temperature history from the 
inside-surface-temperature history. Rearranging equations 
(40) gives 
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m T iA 

Jl 1-0, 

rp _._ T { , 2 -{-d 2 Ti 

12 1 - 0 , 

rp _2 1 <.3'H02?2'1 _ 03^1 

i 3— } a 


y (43) 


rp Ti_ m J rd2Tm~\ J r^2,Tm-‘i J r ■ ■ - + d m Ti 

l-0i J 

The rate of heat flow at any time is_ then determined from 
equation (42). 

APPLICATION 
GENERAL CONSIDERATIONS 

The section on application is devoted to demonstrating 
the solution of two types of problems: In the first type the 
heat-transfer coefficient and adiabatic-wall-temperature 
histories are given and the /wall-temperature solutions are 
obtained. In the second tjqie the temperature history is 
known and the heating-rate history is computed. In each 
example, the problem chosen was a special case, selected so 
that its solution could be and was obtained by an exact 
analytical method. The degree of exactness of the present 
method is demonstrated by comparing each result with the 
solution calculated by exact theory. 

With one exception, radiation was a negligible considera- 
tion in the examples given. Although the method presented 
is well suited to accounting for radiation and includes terms 
for that purpose, radiation was neglected in all cases to make 
possible an exact analytical solution for comparison. 

EQUATIONS FOR HEATED SURFACE 

If the temperature of the heated surface of the wall is 
required, equations (18) and (23), to (25) are used. Whether 
the objective is to compute wall temperature or to compute 
the heat flow from a known wall temperature, the first steps 
involve the determination of the required memory coeffi- 
cients. In order to minimize the labor involved, the recom- 
mended procedure is as follows: First, choose a tentative 
time interval S which seems appropriate to the particular 
problem. (A review of the examples presented herein will 
give an idea of a reasonable value.) Then compute a tenta- 
tive value of the dimensionless diffusion number kS/P. From 
table I or II pick a diffusion number close to the one tenta- 
tively computed. The memory coefficients M and inside 
temperature coefficients 0 given in the table for this diffusion 
number are to be used, and they do not therefore have to 
be computed. An adjustment in the value of 5 is made by 
multiplying it by the ratio of the tabular value of kS/P 
selected to the value of k8/P tentatively computed. Then 
compute from equations (18) and (23) the value of 0 and 
the values of H ; if the radiation is important, B must be 
computed also (eq. (24)). The temperature history of the 
heated-wall surface is then found from equations (25). If 
the heatiDg rate is being determined from a known tempera- 
ture history, the values of H and B are not required. The 
procedure for this case is discussed in the section entitled 
“Example 6.” 


The temperature formulas were derived with the assump- 
tion that the initial wall temperature was zero in order to 
avoid writing T-T 0 numerous times in the formula. The 
simplest way to handle most problems is to subtract the 
amount that the initial wall temperature is above zero from 
both the wall and the adiabatic-wall temperatures. The 
last step in the problem is to add this amount to the solution. 

EQUATIONS FOR INSIDE SURFACE 

If the temperature of the imheated side of the wall is 
required, equations (40) are used. Except for the case of 
the thermally thin wall, all equations were derived for con- 
stant material properties. For small changes in material 
properties with temperature, it appears reasonable to use 
an average value of the properties for the temperature range 
involved. For cases in which material properties vary, it 
seems possible that a more accurate solution might be ob- 
tained by varying the diffusion number or by varying the 
step size to keep the diffusion number constant; however, 
any consideration of such a technique is beyond the scope 
of this paper. 

APPLICATION OF METHOD IN SPECIFIC EXAMPLES 

The following illustrative examples were calculated before 
table I was prepared. Therefore, the values of the coeffi- 
cients M and 0 were computed for the particular walls and 
chosen time intervals S of the examples. All values of M 
and 0 used are listed in table II, which may be considered 
as being supplementary to table I. 

EXAMPLE 1 

Problem for example 1 (a). — A copper wall which is }i 
inch (K< ft) thick is initially at a temperature of 0° F. One 
surface is heated by a boundary layer while the other side 
is insulated. The effective boundary-layer temperature T aw 
is initially 0° F but increases linearly at the rate of 1,000° F 
per second for 10 seconds. The heat-transfer coefficient 
remains constant at A=100 Btu/(hr)(sq ft)(°F). The con- 
ductivity K and diffusivity k of copper are taken as 

K—227 (Btu) (ft) /(hr) (sq ft)(°F) 

&=4.41 sq ft/hr 

Find the temperature history of both wall surfaces. 


Solution for example 1 (a). — The material properties are 
usually givenjin terms of the hour unit. However, since 
fast heating conditions may be more easily understood in 
terms of seconds, time is referred to in seconds and is con- 
verted to hours for use in the equations. For example, if 


5=1 "““seoo 111 


then 


kS (4.41) (24) a 


l 1 


3600 


=0.7056 


By using this dimensionless diffusion number, the values of 
M in column 2 of the following table are obtained from table 
II. The values of T aa are listed in column 4. The value of 
O is given by the equation 



227 

(24) (4.41) 


2.1447 Btu/(sq ft) (°F) 
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TABLE II.— VALUES OF M AND 6 USED IN ILLUSTRATIVE EXAMPLES 


\wi« 

771 

0. 00980 

0. 08820 

0. 17640 

0.35280 

0. 70560 


M 

l 

0. 09187272 

0. 27561815 

0. 38974974 

0. 54911553 

0. 75160198 

2 

-. 01576287 

—.04735497 - 

-. 07101815 

— . 14414264 

-. 35255882 

3 

— . 02656306 

-. 08094552 

-. 13303643 

-. 24244390 

-.33055113 

4 

-. 00947905 

— . 03275298 

-. 06587129 

-. 09447354 

-. 05648228 

5 

00547837 

-. 02295825 

-. 04229223 

-. 03955774 

-. 00990390 

6 

— . 00369637 

01799364 

— . 02736081 

-. 01656451 

-.00173660 

7 

-. 00271361 

— . 01440765 

-. 01770513 

-. 00693626 

-. 00030451 

8 

-. 00210254 

— . 01158046 

-. 01145704 

-. 00290451 

-. 00006339 

9 

-. 00169169 

— . 00931431 

-. 00741389 ' 

-. 00121624 

-. 00000936 

10 

— . 00139980 

-. 00749250 

-. 00479755 

-. 00050929 

-. 00000164 

11 

— . 00118389 

— . 00602714 

— . 00310451 

-. 00021326 

-. 00000029 

12 

— . 00101926 

-. 00484839 

-. 00200894 

-. 00008930 

-. 00000006 

13 

— . 00089082 

-. 00390018 

-. 00129999 

-. 00003739 

-. 00000001 

14 

— . 000788S2 

00313741 

— . 00084123 

— . 00001566 

-. 00000000 

15 

— . 00070668 

— . 00252382 

-. 00054436 

-. 00000656 


10 

— . 00063977 

-. 00203023 

-. 00035226 

-. 00000275 


17 

-. 00058474 

-. 00163317 

-. 00022795 

-. 00000115 


18 

— . 00053905 

-. 00131377 

— . 00014761 

-. 00000048 


19 

-. 00050081 

-. 00105683 

-. 00009545 

-. 00000020 


20 

-. 00046851 

— . 00085014 

-. 00006177 

-. 00000008 



e 

1 

L 00000000 

0. 99289947 

0. 94364543 

0. 80477935 

0. 58038199 

2 

-. 00000012 

-. 09850808 

-. 27773216 

-. 44879472 

-. 47463257 

3 

-. 00002053 

-. 15766415 

-. 23367386 

— . 20690977 

-. 08720676 

4 

-. 00028945 

-. 14162795 

-. 15250956 

-. 08665088 

-. 01629129 

5 

00128807 

-. 11603872 

— . 09871520 

-. 03628441 

-. 00268126 

6 

-. 00314564 

-. 09364233 

-. 06387947 

-. 01519382 

-. 00047015 

7 

-. 00562973 

-. 07537042 

-. 04133602 

-. 00636230 

-. 00008244 

8 

-. 00802789 

-. 06063594 

— . 02674906 

00266417 

-. 00001440 

9 

-. 01036556 

-. 04877804 

-. 01730940 

-. 00111560 

-. 00000253 

10 

-. 01237062 

-. 03923850 

-. 01120097 

-. 00046715 

-. 00000044 

11 

-. 01402862 

-. 03156463 

-. 00724818 

-. 00019562 

-. 00000008 

12 

-. 01533940 

-. 02539137 

-. 00469032 

-. 00008191 

-. 00000001 

13 

-. 01633809 

-. 02042651 

-. 00303512 

-. 00003430 

00000000 

14 

-. 01706855 

-. 01643084 

— . 00196-104 

-. 00001436 


15 

— . 01757475 

-. 01321741 

-. 00127093 

— . 00000601 


16 

— . 01789689 

-. 01063245 

— . 00082243 

— . 00000252 


17 

-. 01806993 

-. 00855303 

-. 00053219 

-. 00000105 


18 

-. 01812340 

-. 00688029 

-. 00034438 

-. 00000044 


19 

-. 01808182 

-. 00553469 . 

-. 00022285 

-. 00000018 


20 

-. 01796529 

-. 00445226 

— . 00014421 

-. 00000008 



The value of /t=100 Btu/(hr)(sq ft)(°F) is converted to 
.5=0.00800. The use of columns 2 and 4 in equations (25) 
gives the heated-wall temperature in column 5. Using the 
inside-surface-temperature formulas, equations (40), and 
the values of 6 in column 6 gives the values of T t in column 7. 


1 

2 

3 

4 

5 

6 

7 

Term 

M 

Time, 

T aa , 

T, 

e, 

T,, 

number 


sec 

°F 

°F 

°F 


0 


0 

0 

0 


0 

1 

0. 75160 

1 

1, 000 

11 

0. 580382 

4 

2 

-. 35256 

2 

2, 000 

36 

— . 474632 

21 

3 

-. 33055 

3 

3, 000 

74 

— . 087207 

50 

4 

-. 05648 

4 

4,000 

124 

-. 015291 

91 

5 

-. 00990 

5 

5, 000 

186 

-. 002681 

144 

6 

-. 00174 

6 

6, 000 

260 

— . 000470 

210 

7 

-. 00030 

7 

7, 000 

346 

-. 000082 

288 

8 

-. 00005 

8 

8,000 

443 

— . 000015 

377 

9 

-. 00001 

9 

9, 000 

552 

-. 000002 

478 

10 

-. 00000 

10 

10, 000 

073 

-. 000000 

591 


The wall-surf ace-temperature curves of T and T t are 
shown in figure 1 (a). For comparison, the results calcu- 
lated by the theoretically exact formula (eq. (B10) of ap- 
pendix B) are shown. This formula is 


72 eo 

r= 6 f- 26 ^s 

1C n- 1 


cos(^ a ^){l~e pan ) 

a "( 1 +! K + fl a ') COSa * 


(44) 


where b represents the slope of the adiabatic-wall-tempera- 
ture curve and a n represents the positive roots of the auxil- 
iary equation 

a, tan a„=^ (44a) 

Setting 6=3,600,000 °F/hr and setting x=l and *=0 in 
equation (44) result in the values for outside and inside 
temperatures plotted as circles and squares in figure 1. 
The comparison shows that accurate results are obtained 
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(a) H-inch-thick wall. 

Figure 1. — Example 1. Temperatures of copper-wall surfaces. 
Adiabatlc-wall temperature varies linearly from 0° to 10,000° in 
10 seconds; ft = 100 Btu/(hr)(sq ft)(°F). 


by the present method, the maximum difference between 
methods being about 1°. 

Problem for example 1 (b). — The conditions for example 
1 (b) are the same as those for example 1 (a), except that 
the copper wall is 3 inches thick, or l=){ foot. 

Solution for example 1 (b). — Since this wall is so thick, 
more highly transient conditions prevail throughout the 
heating period. A value of 5 of K second or 1/7200 horn 
was therefore used. With the same procedure as used be- 
fore, the results of the present method are given as con- 
tinuous lines in figure 1 (b), while the results from equation 
(44) are given by the symbols. The agreement is evident. 

EXAMPLE 2 

Problem, — Example 2 illustrates the principal advantage 
of the present method; that is, its capability of handling 
arbitrary variations of T aa and h. Neither is it necessary 
to know a mathematical formulation for these variations. 

A K-inch copper wall which is initially at a temperature 



Figure 1. — Concluded. 


of zero is aerodynamically heated on one side and insulated 
on the other. The time histories of T aw and h are given at 
^-second intervals in the following table: 


Time, 

sec 

T„, 

°F 

h, Btu/(hr) 
(sq ft) ( °F) 

Time, 

sec 

I'atey 

°F 

h, Btu/(bx) 
(sq ft) ( °F) 

0 

0 

36 

5. 5 

5, 391 

68. 4 

. 5 

1, 365 

41. 4 

a 0 

5,366 

66. 6 

1.0 

2, 485 

45. 0 

a 5 

5,255 

63. 9 

1. 5 

3, 388 

48. 6 

7. 0 

5, 107 

60.0 

2. 0 

4, 094 

52. 2 

7. 5 

4, 831 

55. 8 

2. 5 

4, 620 

55.8 

a 0 

4,335 

52. 2 

3. 0 

4; 932 

60. 0 

& 5 

3, 654 

48. 6 

3. 5 

5, 119 

63. 9 

9. 0 

2, 769 

45. 0 

.4. 0 

5,263 

66. 6 

9. 5 

1, 658 

41. 4 

4. 5 

5, 345 

68. 4 

10. 0 

297. 5 

3a 0 

5.0 

5, 387 

69. 0 





The temperature history of both wall surfaces is to be 
found. 

Solution, — The heating conditions are severe and con- 
tinuously transient, with the boundary-layer temperature 
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rising and falling over 5,000° in 10 seconds. A computing 
interval 5 smaller than that used in the first problem is 
therefore used. Let S=K second= 1/7200 hour. The re- 
sulting wall -surface-temperature curves are drawn in figure 
2. The circles and squares give the surface temperatures 
computed by a theoretically exact procedure. Comparison 
shows the present method to be accurate. 

EXAMPLE 3 

Problem. — Example 3 is the same as example 2 except 
that the wall is 3 inches thick and the effective boundary-layer 
temperatures are slightly different. The initial value of 
T aw is 0. The subsequent values of T am are given at 0.5- 
second intervals by the following time series: T aw = 1,365, 
2,484, 3,386, 4,088^ 4,609, 4,915, 5,094, 5,227, 5,296, 5,325, 
-5,315, 5,265, 5,149, 4,98.6, 4,694, 4,184, 3,489, 2,591, 1,471, 
100. Find the temperature history of both wall surfaces. 

Solution (a) (thick-wall solution). — The value of S was 
taken as ]i second. By following the procedure previously 



O 2 4 6 8 10 

Tima, sec 

Piguke 2. — Example 2. Temperatures of )4-inck copper wall heated 
according to assigned history of h and T aw - 


presented, the wall-surface-temperature curves shown in 
figure 3 were obtained. Again the symbols give the results 
of exact theory. 

It should be noted that the “exact solution” for examples 
2 and 3 is not actually an alternate method of solution for 
any practical problem but gives a solution to the particular 
problems only. The “solution” was obtained by working 
in reverse; that is, a heat flow was assumed and the corre- 
sponding boundary-layer characteristics were computed. 
A truly alternate method of solution is now considered, 
however. 

Solution (b) (infinitely thick-wall solution).— Since the 
thermal lag of a 3-inch copper wall is so great when subjected 
to the rapid heating specified by this problem, it appears 
reasonable to obtain the heated-surface temperature by 
assuming that the^wall is infinitely thick. The memory 
coefficients are the same for all infinitely thick walls and are 
given by equations (31). The same temperature formulas, 
the same values of T am , and the same values of h are used 
as before, but the values of are given by equation (29). 
The results calculated by this method are listed along with 



Figdbe 3. — Example 3. Temperatures of 3-inch copper wall heated 
according to assigned history of h and T aw . 
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those from the thick- wall solution (a). 



T, °F, calculated by — 

Time, 



sec 

Thick-wall 

Infinitely 


solution 

thick-wall 



solution 

i 

12. 8 

12. 8 

2 

35. 3 

36. 3 

3 

01.3 

6L 3 

4 

87. 0 

87. 0 

6 

109. 1 

109. 1 

0 

124. 8 

124. 8 

7 

132. 3 

132. 3 

8 

130. 0 

130. 0 

9 

119. 5 

119. 4 

10 

100.0 

99. 8 


The consistency of the alternate methods for a ther m all y 
thick wall is evident. The reason for the close agreement, 
may be found in figure 3, which shows that the unheated 
surface of the 3-inch wall rose to only 14° F. 

EXAMPLE i 

Problem. — The most severe test of the present method 
would occur if there were a large instantaneous increase of 
T aa . While this condition could hardly happen in flight, it 
might happen if a research model were suddenly immersed 
in a high-stagnation-temperature jet. Let a #-inch copper 
wall, initially at a temperature of zero, be instantly subjected 
to an effective boundary-layer temperature of 5,000° F on 
one surface while no heat transfer occurs on the other sur- 
face. The heat-transfer coefficient is 100 Btu/(hr)(sq ft) 
(°F). Solve for the temperature history of both wail sur- 
faces for 10 seconds. 

Solution. — In this case not only is there a very high 
transient-temperature condition initially but the instanta- 
neous increase in T aa does not lend itself to approximation 
by the unit triangle. The simplest procedure is to take 
small steps for the first few seconds to minimize the errors 
introduced. In order to help circumvent the difficulty of 
calculation, an excellent method of approximating the wall 
surface temperature for the first or first few small steps is to 
use the following formula from page 109 of reference 5, 
which gives the temperature on the surface of an infinitely 
thick wall for a constant flow of heat at the surface: 



Time, sec 


Figure 4.— Example 4. Temperatures of J4-inch copper wall 
after application of 5,000 °F jump in gas temperature. ft=100 
Btu/(hr) (sq ft) (°F) . 


define the highly transient part of the curve, then JLsecond 
steps for the remaining 9 seconds. 

A theoretically exact solution to this problem was obtained 
by equation (B7) of appendix B. The results of applying 
this equation are given by the symbols in figure 4, which 
shows that agreement was obtained. 


EXAMPLE 5 



The values of 5 used were S=% second for 2 seconds, then 
5=% second for 8 seconds. Since the use of an equation 
based on an infinite wall is permissible for a K-inch copper 
wall for at least / second, the values of T were computed by 
equation (46) for the first two ^-second steps, then by the 
usual equations. The results are presented in figure 4. 
Since the inside temperature T t depends on the outside 
temperature T and not directly on T aa , there is no particular 
difficulty of approximation in obtaining T t . Accordingly, 
in obtaining T it ’(-second steps were taken for 1 second to 


Problem. — A Ke-inch Inconel wall is heated by a high- 
temperature jet. If T aa — 5,000° F and h=50 Btu/(hr)(sq 
ft)(°F), determine the skin- temperature history for 15 sec- 
onds. Neglect radiation. 

Solution. — The heat capacity of a fis-inch Inconel wall was 
assumed to be <?= 0.3229 Btu/(sq ft) (°F) . The example was 
worked three times with values of S of 1, 2, and 5 seconds to 
show how sensitive the ‘Thin wall” formula is to the time 
interval. Substituting the given constants into equation 
(34) yielded the results given by the symbols in figure 5. 
For this example the exact theory is shown by the solid line. 
Evidently, large time intervals are permissible with this 
formula. 
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Time.sec 

Figure 5. — Example 5. Temperatures ofDlo-inch Inconel wall. 

7aw=5,000 °F ; A=50 Btu/(hr)(sq ft)(°F). 

EXAMPLE 6 

Problem for example 6 (a). — The temperature history of 
the heated surface of a %-inch copper will initially at zero 
temperature is given by the following time series in which 
the temperatures are for ^-second intervals: T— 4.8, 13.2, 
25.2, 40.2, 58.2, 77.9, 99.4, 122.6, 146.3, 169.8, 193.0, 214.8, 
235.3, 253.2, 268.8, 281.5, 290.8, 296.7, 298.0,' 297.5. The 
inside surface is insulated. Determine the history of heat 
flow into the heated surface from the given surface-tempera- 
ture history. 

Solution for example 6 (a). — If the time interval used is 
sufficiently small, the average rate of heat flow over the 
interval is a good approximation to the rate of heat flow 
at the center of the interval. Equation (42), which gives 
the average rate of heat flow over the interval ending at 
t=m5, may be used. For a value of 5 of % second, the 
values of M are given in table II . Substituting in equation 
(42) gives the rate of heat flow plotted as circles in figure 6. 
The solid curve gives the theoretically exact instantaneous 
rate of heat flow for comparison. The results from equation 
(42) are seen to be precise. The system yielding instan- 
taneous heat flow, mentioned previously, would seem to he 
a natural one for the present problem; however, the results 
obtained by that system were found to be inferior to those 
presented. 

Problem for example 6 (b). — The corresponding inside- 
surface-temperature history of the same wall is given by 



Figure 6. — Example 6. Rate of heat flow into heated (outer) sur- 
face of }4-inch copper wall computed from temperature history of 
outer surface and from temperature history of inner surface. 


the following time series in which the temperatures arc for 
^-second intervals: T { = 0.8, 4.4, 11.7, 22.4, 36.0, 52,6, 

71.6, 92.3, 114.2, 137.4, 160.6, 183.7, 205.4, 226.4, 245.3, 

261.6, 275.2, 285.6, 292.4, 296.3. Determine the history 
o'f heat flow into the outer surface by using only the given 
inside temperatures. 

Solution for example 6 (b). — Equations (43) may be used 
to determine the outside-surface-temperature histoiy from 
the inside-temperature history. Then, the rate of heat 
flow at any time is determined as in the solution for example 
6 (a). The factor 1/(1— 0 X ) in equations (43) may bo 
thought of as a magnification factor. Large values of this 
factor tend to cause an instability in the computed tempera- 
tures. In this example, if 5 is taken as K second or 1 /7200 hour, 
the value of kS/l* is 0.3528, 8i is 0.80478, and the magnifica- 
tion factor is 5.1. On substitution in equations (43), an 
oscillation of period 26 builds up in the values of T, veiy 
slowly at first and to either side of the correct answer as a 
mean value, but in a divergent manner so that by 4 seconds 
the amplitude is 6°. For a value of 5 of 1 second, the value 
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of kSfP is 0.7056, 6i is 0.58038, and the magnification factor 
is 2.4. In this instance, an oscillation in T of period 25 and 
maximum amplitude of 1.5° occurred. Substituting these 
values of T, without fairing, into equation (42) gives the 
results shown by the square symbols in figure 6. If a larger 
value of 5 were used, the oscillation would be damped out 
but the accuracy would suffer because of a lack of definition 
of the rapidly Varying heating rate. The particular case 
demonstrated is therefore a marginal one for the determina- 
tion of heat flow from the temperatures of the inside surface. 
The instability is found to disappear for thermally thinner 
walls and conversely to increase rapidly for thicker walls. 


CONCLUDING REMARKS 

Formulas to facilitate the determination of the transient 
surface temperatures of thick walla from an arbitrary 
variation of adiabatic-wall temperature and heat-transfer 
•coefficient have been developed. Formulas to facilitate the 
determination of heat flow from an arbitrary variation of 
wall surface temperature were also obtained. The numerical 
applications given demonstrate a high degree of accuracy 
for the present method. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 18, 1957. 



APPENDIX A 

REVIEW OF SUPERPOSITION AND TIME SERIES 


TIME SERIES 

A time series may be defined as a series of numbers or 
quantities which represent, the values of a function of time 
at successive equal intervals of time. According to the 
notation of Tustin (ref. 3), each quantity is separated from 
the others by a comma since the values corresponding to 
different times are not added together. The quantity at 
zero time is zero. The first quantity recorded represents 
the value of the function at the end of the first time interval. 
The time interval used is arbitrary and its size is selected 
to obtain the accuracy required in the solution of a specific 
problem. The symbol for the time interval is 5. Thus, the 
series y=y i} y 3 , y 3 , . . . y m represents the values of the 
function y at the times 5, 28, 3 8, .. . m8. 

THE UNIT TRIANGLE 

A unit triangle is an isosceles triangle which has an altitude 
of unit magnitude and a base of 28, or two time intervals. 
Since 5 is an arbitrary time interval, the unit triangle is 
accordingly arbitrary. A plot of a unit triangle centered 
at t—8 is shown in sketch 1, where y represents magnitude 



Sketch l. 


or altitude. The slopes of the sides of the triangle depend 
on the value of 8 and are equal to ± 1/5. Three lines may be 
superimposed to represent the unit triangle. The equations 
of the three lines which may be added to represent the unit 
triangle are y=(l/S)t, y= — (2/8)(t—8), and y=(l/5)(£— 25). 
Although the triangle terminates at t=28, the values of £ in 
the equations for y can go on to infinity since the values of y 
add to 0 beyond £=25. 
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THE FUNCTION CORRESPONDING TO A EEFEBENCE LINE 

Let F(x,t) be the solution to a boundaiy-value problem 
specified by a linear partial differential equation and the 
linear boundary condition y— (l/5)£, where y is the value of 
F or one of its derivatives or integrals at some fixed value of 
x. Because of the linearity of the problem, the magnitude 
of the solution is directly proportional to the magnitude of 
the slope 1/5 of the boundary condition. For example, if 
y= (2/5)£, the corresponding function representing the solu- 
tion is 2 F(x,t). Again, if y— — (2/8)t, the corresponding 
function is —2 F(x,t). The slope of the line y= (l/5)£ can 
thus be used as a reference magnitude. This slope is the 
same as that of the left side of the unit triangle. 

The value of £ in F(x,t) is always identical with the value 
of £ in the boundary condition y=(l/8)t. Thus, if the origin 
is shifted so that t/=(l/5)(£— 25), then the corresponding 
function is F(x,t—28). Particular solutions of a linear 
differential equation can always be added in linear combi- 
nations to satisfy more general boundary conditions. If a 
and b are constants, and aF(x,t) corresponds to the bound- 
ary condition y l ={al8)t, and bF(x,t—8 ) corresponds to 
Vi— (bl5) (t—8), then the function corresponding to the sum of 
the two lines y=yi+y 3 is F—aF(x,t)+bF(x,t—8). Let an 
additional property of F, as well as of y, be that it assumes 
the value 0 for any time less than 0. The range of time of 
interest is therefore from 0 to °° . 

THE FUNCTION CORBESPONDING TO A TRIANGLE 

Now consider a function F&(t) dependent on the lines of 
the unit triangle in sketch 1 for its value just as F(t) is 
related to the line y=(l/5)£. In consideration of the three 
lines by which the unit triangle may be replaced, y—(\l8)t, 
y= — (2/8) (t—8), and y= (1/5) (£—25), the three correspond- 
ing solutions or functions of time are F(t), —2 F(t—8), and 
F(t—28). Because of the additive nature of solutions, the 
solution corresponding to the complete triangle may be 
defined as the sum of the solutions for the lines which com- 
pose it: 

F*(t)=F(t)-2F(t-8)+F(t-28) (Al) 

In equation (Al) each term has a value of zero for negative 
values of its argument. 

BEPBESENTATION OF A GENERAL CURVE BY TRIANGLES 

The curve A-B in sketch 2 (a) is any arbitrary continuous 
function In the plane y,t which can be faired through its 
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ordinate values 0 , y lt y 2 , and y 3 . The curve may be well 
approximated by a series of chords such as a, b, and c if a 
suitable spacing 5 is used. 




Sketch 2. 

Draw the lines d, e, f, g, and h of sketch 2 (b) to form the 
isosceles triangles with sides a and d, e and f, and g and h. 
Since the sum of two straight lines is a straight line, it is 
clear that, if the letters which designated the lines are now 
used to designate the equation of the line, then d+e=b and 
f-f-g=c. Of course, the line a is the first half of the first 
triangle as well as the first chord of the curve A-B. In 
designating the ordinates y u y 2 , y 3 , . . . y m , the altitudes of 
three triangles whose sides add up to the chords of the curve 
A-B are simultaneously designated. The ordinate series 
Hu Vi> thy ■ ■ • Vvx> with spacing 5, is a time-series approxi- 
mation of the curve A-B. In this case, as well as elsewhere 
in this report, each ordinate of a series is understood to 
represent the altitude of an isosceles triangle with a base 
width of 25. 


AKBITBABY CONTBOL FUNCTIONS 

Suppose that the curve A-B, or y(t), is to serve as an 
arbitrary reference or control function for a corresponding 
function F v (t). It is desired to determine F v (t) in terms of 
functions corresponding to unit triangles, such as F&(l), 
given by equation (Al). The function F&(t) is a function 
corresponding to a unit triangle centered at t=8. If F\(f) is 
multiplied by y u there results the function corresponding to 
the triangle with sides a and d of sketch 2 (b). By super- 
position, the function corresponding to the three triangles 
centered at 5, 25, and 35 is 

F y (t) =y x F& (t) + ijiF^ (t— 5) 4 -yJF i (t — 25) (A 2 ) 

Then, since the three triangles add up to the chords of the 
curve A-B, this is the function, in ordinary algebraic form, 
corresponding exactly to the chords of the general reference 
curve y=yi, y%, y 3 , . . . y m . In order to put equation (A 2 ) 
in time-series form, let Fiy(t)=a 1) ch, a 3 , . . . a m , with spacing 
5. Substituting in (A 2 ) and placing terms for the same time 
in columns and adding yield 

ViFty(t) =yia u t/iOs, z/i®s» . . . 

ViF^it—S) = 0, y& 1 , y 2 da, . . . 

ysFiit— 25)= 0, 0, y& u . . . 

Fy(t) =yio-i, (yich+y&i), (z/i^+z/^+z/aUi), • • . (A3) 

The result shown by equation (A3) is obviously that which 
is obtained by formal algebraic multiplication of a u a-i,a 3 , . . . 
by z/i, Vi, y 3 , ■ ■ ■ as follows: 


F&(t)= Oj, 02, 03, 

y(t)= vu y-i, y>, 

ViO-ii Z/i®2, ViChy 

y&Xy Vtfht 

ys(h, 


Fyity—yyd!, (yiOi+yufli), (z/iO^ +Z/ 2 O 2 + Z/yz 1 ) , • . . (A 4 ) 
Therefore, 

Fy(t),=[F^m[y(f),] (AS) 

where the subscript 8 denotes time-series form. 

Equation (A5) states a simple theorem which was first 
given in reference 3. It was developed with. y(t) having a 
value of 0 at f =0 and F&(t) also having a value of 0 at t= 0 . 
If either or both of the series had values other than 0 at 
t= 0 , multiplication as in equation (A5) would not be suffi- 
cient to obtain F v . 



APPENDIX B 

SUMMARY OF ANALYTICAL TEMPERATURE FORMULAS FOR THICK WALLS 


CONSTANT FLUID TEMPERATURES 

Consider an infinite wall of thickness l and an initial tem- 
perature of zero. (See sketch 3.) Let the wall be suddenly 



contacted at the face x—l by a fluid of temperature T aa while 
the face at x=0 remains exposed to a fluid of zero tempera- 
ture. Let the heat-transfer coefficient at the face x=l have 
the value h and the heat-transfer coefficient at the face x=0 
have the value rih . Let the wall have uniform physical prop- 
erties which are invariant with time. 


The flow of beat within the wall is governed by Fourier's 
equation for transient heat flow, which states that the rate 
of increase of temperature is proportional to the rate of 
change of temperature gradient: 

l 

bT b v T 
bt~ k bx* 

The constant of proportionality, called the diffusion coeffi- 
cient, is equal to the ratio of the conductivity of the material 
to the heat capacity as represented by the product of specific 
heat and weight density as follows : 




The boundary equation for the wall at x—l is obtained by 
equating the heat transfer between the fluid at temperature 
T au> and the wall at temperature T to the rate of heat transfer 
in the wall at x=l: 


h[T at0 —T x .i)=K (B2) 

A similar boundary equation is written for x—0: 

D„. (B3) 

The initial condition specifies that the initial temperature of 
the wall is zero : 


(T)<- o=0 (B4) 


'I'he simultaneous solution of partial differential equations (Bl) to (B4) by operational methods yields the infinite sories: 


T=T ai 


, , hx 
1+71 K 

l+’7+’7 ^ 


-2 T att S 



cos (f«,)+igsm (£«,)_ 

e * 

n k . hi 
L hl an v K 

-3(1 + 77) 

cos a *+[0^+ 1 +7 i 

\ 9 v hr 

sin a* 


(B5) 


The parameters a* are angles which are the positive roots of 
the equation 

tan ot„= a jz +v) u (B5a) 

an hl~ v K 

Two dimensionless nmnbers of physical significance are 
present in equation (B5) : the diffusion number kt/l 1 , and the 
conductance or Nusselt. number hl/K. The first term of 
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equation (B5) gives the steady-state solution or equilibrium 
condition. 

' SPECIAL CASE ,-0 

The special case of a plate that is heated by convection 
on one face while the heat transfer on the other face is so 
small as to be negligible is important. Setting the heat- 
transfer coefficient at the unheated face equal to zero cor- 
responds to making the assumption that the plate is per- 
fectly insulated at that face. If the plate is perfectly in- 



A METHOD OF COMPUTING THE TRANSIENT TEMPERATURE OF THICK "WALLS 


895 


sulated at the face x= 0 , substitute 17=0 in equations 
(B6) and (B6a) and utilize the following relation to eliminate 
sin a B , 

hi 


an tan a a = 


K 


(B6) 


where a n represents the positive roots of the equation. 
These substitutions yield the following expression for tem- 
perature : 


T=T ao - 2 T au T, 


cos (j 


Cx \ -4 a.* 


r«) 


\ l+ K + U an ) 003 a * 


(B7) 


Equations (B6) and (B7) may be found on page 100 of ref- 
erence 1. 

SPECIAL CASE ij-0 AND h= « 

The special case in which ij=0 and h= corresponds 
physically to perfect insulation on one face and a known 
initial temperature on the other face. This case is developed 
by setting h= <» in equations (B6) and (B7). Equation 
(B6) becomes 

a* tan a„= <*> 
tan a n = a> 

or, since a* assumes the sequence -g Vjr » ~2’ an< ^ 80 f orfc h, 


a 0 =(2n— 1) | 
Equation (B7) reduces to 

r=T »-; r \liSr c 09 [ ( 2 '^ 1 ) 5 f] 




(B8) 


Equation (B8) may be found in standard references (ref. 5, 
page 196, problem 6, for example). In this equation, T a 
is the stop in wall-surface temperature used in place of 
Tau,. 

VARIABLE FLUID TEMPERATURE 

ARBITRARY VARIATION OF T ca 

Equation (B7) is an exact solution of a thick-wall hound- 
ary-layer heating problem which is suitable for checking 
the accuracy of the present method of computing wail 
temperatures. Another method, which is more general, 
can be obtained by letting the adiabatic-wall temperature 
vary in a known manner, T aw =F(t). Let V(x,t) be the 
variation of wall temperature due to a 1° step in adiabatic- 
wall temperature. If the initial wall temperature is zero, the 
wall temperature is given by Duhamel’s formula as 


T(x, t)= j‘ F(t)V{x, t-T)di 


m 


where F is the derivative of F with respect to t and r is a 
dummy time variable. 


LINEAR VARIATION OF T 


In the linear variation of T aa , T aa =bt, or, in the notation 
of Duhamel’s formula, F=bt and F=b. Substituting T 
from equation (B7), with T aa = 1, into equation (B9) and 
performing the integration give 


0172 « cos 
T(x,t)=bt-^J1 . 

rC n - 1 




(BIO) 


COS «n 

Equation (B6) is used to obtain a„. 

LIN E AR VARIATION OF WALL SURFACE TEMPERATURE 

One of the fundamental equations on which is based the 
time-series development of surface-temperature equations is 
one expressing the transient wall temperature due to an in- 
crease of 1° in temperature of one wall surface during each 
interval of time 5. The assumptions are made that the 
initial temperature is zero and that one wall is insulated. 
From these considerations, F in Duhamel’s formula is 
F=t/S and F=l/S. Substituting this derivative and T from 
equation (B8), with T s — 1, into equation (B9) and perform- 
ing the integration give the following equation for the* tem- 
perature at any plane due to a surface temperature rise of 1° 
per unit of time 5 : 


d 7T 3 kS^i 


(-D 


n— 1 


( 2 n — l) 3 V 1 


x 3 k 

-Vn-l)iL-±t 
„ I p 


(Bll) 


The heat flow at any point within the wall due to 1° rise of 
wall surface temperature in the time 5 is obtained by multi- 
plying the conductivity by the temperature gradient. Dif- 
ferentiating equation (Bll) with respect to x and multiplying 
by K yield 


8 ik ^ 8111 [c^- 1 ) If] ( 

2 IT* *5 £={ (2» — l) 2 V 1 


X J i 

-(2n-l ).— -t 

» I P 


In order to el imin ate k in the coefficient, substitute for k its 
definition k=K/pc. Then lK/k 5 =pcl/ 5 . The product pci 
is the heat capacity of the wall per unit area per degree and 
is represented by the symbol 6 . In order to obtain the heat 
flow at the heated surface, where the heat balance is to be 
made, let x=l. Hence, the instantaneous heat flow due to a 
uniformly increasing surface temperature of 1° per unit of 
time 5 is simplified to 


_ 80 ^ 1 -e 
2 7r>5 tk ( 2 n- 1) ! 


(B12) 
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In order to find the average flow of heat over the interval 
5, equation (B12) is integrated with respect to time between 
the limits (f— 5) and t. This integration gives the total heat 
flow through the surface dining the interval. On dividing by 
5, the average rate of flow for the interval is obtained. The 
result is 


- 

2- tt* 5 § 


1 


t 

— (2n— 1) J — — I 


* *•* k 

-C2»-l)i— (<-*) 


(2n- l) a 


j~S(2n-iy 


^S(2 n -iy 


(B13) 


The summation X) 75 - — yi is a constant with the value 7 ^/ 8 . 

n-l (All lj 

In working with time series only integral increments of time 
(5, 25, 35, . . . 77 / 5 ) are used. Inequation (B13), therefore, 
t may have any value mb, where m is an integer. For con- 
venience, and to systematize results, the following identities 
are defined: 


4 1* ® i 

i § (2n-iy 


-(m-l) 

a - e 

m - 1_ ^5jfe J 1 ( 2n-l )« 

x* £ 

-TO(2n-l)2 — — a 

= _4P_ ~ £ 4 ' 

(27/. — l) 4 


Equation (B13) then becomes 

5 _ _ 


•7T 3 6 \8 


-dm— 1^ 


(B14) 


(BIS) 


for the time interval ending at t=mS. Equation (BIS) is 
used in forming the heat balance at the heated wall surface 
in the derivation of wall temperature formulas. 
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ANALYTICAL TEMPERATURE DISTRIBUTION 


Equations (40) are the formulas most convenient for 
determining the temperature of the inside surface or of any 
plane within the surface. For the determination of the 
temperature of planes within the wall, equations (36a) are 
ordinarily evaluated numerically to obtain 6 for a given ktjP 
and x/l. However, it may be necessary to obtain the tem- 
perature distribution through the wall at some instant of 


time due to an arbitrary surface-temperature history in a 
form that is analytical with respect to x and that can be 
differentiated or integrated with respect to x (for example, 
in the derivation of a general formula for the thermal stress 
distribution or maximum stress). This temperature dis- 
tribution is obtained by substituting equations (36a) into 
equations (39), substituting equations (39) into equations 
(40), and collecting terms. The results are 
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)e 
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X* JW 
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( — T m -i-\-2T„-3 — T m 




+ (2T 1 -T 2 )e 


-( CT -l) (2n-l)»£^ 


4 P 


T,e 






Note that T u u T it 2 , T t , S) and Ti, m can be computed inde- labor, however, because all terms in the summation must be 
pendently. The use of these equations involves considerable summed in unison; 


REFERENCES 


1. Cnralaw, H. S., and Jaeger, J. C.: Conduction of Heat in Solids. 

The Clarendon Press (Oxford), 1947. 

2. Schmidt, Ernst: Thermodynamics. The Clarendon Press (Oxford), 

1949, pp. 428-436. 

3. Tustin, A.: A Method of Analysing the Behavior of Linear Systems 

in Terms of Time Series. Jour. Inst. Elec. Eng. (British), 


vol. 94, pt. n-A, no. 1, 1947, pp. 130-142. 

4. Bubb, Frank W.: Cardinal Spectrum Analysis of Physical Systems. 
AF Tech. Rep. No. 6300, “Wright Air Dev. Center, U. S. Air 
Force, Mar. 1961. 

6. Churchill, Ruel V.: Modern Operational Mathematics in Engineer- 
ing. McGraw-Hill Book Co., Inc., 1944. 


897 




